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From a Dalitz plot analysis oi B —^ Kn^^p' decays, we find a signal for Z(4430)^ -K^ip' 
with a mass M = (4443li2il3) MeV/c^ width V = (lOTtlgtge) MeV, product branching fraction 
B{B° R- Z{U30) + ) X B{Z{U30)+ tv+iIj') = (3.2+J;^+f j) x 10"^ and significance of 6.4 a that 
agrees with previous Belle measurements based on the same data sample. In addition, we determine 
the branching fraction B{B° -* K*{892)°ij') = {5.52+°o ll+°o fg) x 10"^ and the fraction of JC*(892)° 
mesons that are longitudinally polarized /i = (44.8l2'°l5'3)%- These results are obtained from a 
605 fb~^ data sample that contains 657 million BB pairs collected near the T{4S) resonance with 
the Belle detector at the KEKB asymmetric energy e^e~ collider. 

PACS numbers; 14.40.Gx, 12.39.Mk, 13.25.Hw 



INTRODUCTION 

In a paper on the B meson decay process B 
KiT'^ij}' [H, the Belle Collaboration reported the ob- 
servation of a distinct and relatively narrow peak in the 
■K^ijj' mass spectrum near M(7r+V'') — 4430MeV/c^. 
The analysis was performed by excluding the events in 
the M{Ktt+) regions of the A'*(892) and K*{U?,{)), and 
fitting the one-dimensional M(7r+-0') distribution. The 
fit gave a resonance mass and width of M = (4433 ± 
4 ± 2) MeV/c^ and V = (45ti^t?i^) MeV, where the first 
uncertainty is statistical and the second is systematic; 
the significance of the resonance was 6.5cr. If this peak, 
called the Z(4430)+, is interpreted as a meson state, then 
it must have an exotic structure; its minimal quark con- 
tent is \ccud). 

The Z(4430)+ observation motivated a subsequent 
Belle study of the process — > K~Tr^Xci, where a dou- 
bly peaked structure was observed in the tt+^cI invariant 
mass distribution Q . In this channel, the observed struc- 
ture is rather wide, therefore a full Dalitz plot analysis 
was used in order to establish that the observed peaks 
could be unambiguously associated with dynamics in the 
TT+Xci channel. If these peaks, called the Z(4040) and 
^(4240), are attributed to meson states, a minimal four- 
quark substructure is required. 

A recently reported study of B ^ Kn^^p' decays by 
the BaBar Collaboration [J] did not find a significant 
signal for Z(4430)+ 7r+-0'; the reported significance 
is at the 1.9 a — 3.1 a level. The BaBar sample of i? ^ 
Kn^ip' decays is about 85% the size of the corresponding 



Belle data sample. 

In this paper we present a reanalysis of the Belle B 
Kn'^ip' data sample using a Dalitz plot formalism. 

The Belle detector Q is a large-solid-angle magnetic 
spectrometer that operates at the KEKB asymmetric- 
energy e~^e~ collider A data sample corresponding 
to an integrated luminosity of 605 fb~^ collected at the 
T(4S') resonance and containing 657 million BB pairs is 
used. A GEANT-based Monte Carlo (MC) simulation 
is used to model the response of the detector. 



EVENT SELECTION 

We select events of the type i?° K~tt^iI'' and 
_B+ KgTT^ip', where the ip' decays either to or 
TT+TT-J/iP with J/iP £+£- = e or /x), K°_-^ tt+tt-. 
We use the same selection criteria as in Ref. [21 • In par- 
ticular, we identify B mesons using the beam-energy con- 
strained mass Mbc = \/ £'bcam ^ P% energy dif- 
ference AE ~ -Ebcam ~ Eb, where £'boam is the center- 
of-mass (cm.) beam energy, pB is the vector sum of 
the cm. momenta of the B meson decay products and 
Eb is their cm. energy sum. We select events with 
A/bc — ttibI < 7.1MeV/c^ (m^ is the world-average B- 
mcson mass 0) and \AE\ < 34 MeV, which are both 
±2.5(7 windows around the peak values. To model com- 
binatorial backgrounds, we use events that are in the 
A/bc signal region and the AE sidebands defined as 
\AE ± 70 MeV I < 34 MeV. To improve the definition 
of the Dalitz plot boundaries for both signal and side- 
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band events, we perform a mass-constrained fit to the 
B candidates from both regions. Simulations of the two 
tp' decay modes indicate that the experimental resolution 
for M{Tr+ip') is cr = 2.5MeV/c2 for both modes. 

DALITZ PLOT DISTRIBUTION 

We sum the Dalitz distributions for _B" K^n^ip' 
and KgW'^il'' candidates. Due to the mass differ- 

ence between K~ and Kg the corresponding Dalitz plots 
have slightly different boundaries. We find that this has a 
negligible effect on the results of the Dalitz analysis. The 
Dalitz plot for the AE signal region is shown in Fig. [T] 
Here vertical bands corresponding to the K*{892) and 
the K* (14:30) are evident. The horizontal cluster of en- 
tries in the vicinity of A'P{7r+ilj') ~ 20GeV^/c'* consti- 
tutes the Z(4430)+ signal reported in Ref. [2|. 
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FIG. 1: The B Kn+i)' Dalitz plot for the AE sig- 
nal region. The solid (dashed) lines delimit the five vertical 
(three horizontal) slices that are used to present the fit re- 
sults in subsequent figures. The coordinates of the vertical 
lines are M''{Ktv+) = (0.796)^ GeVVc*, (0.996)^ GeVVc^ 
(1.332)2 GeVVc* and (1.532)^ GeVVc"; the coordinates of 
the horizontal lines are M'^{tv'^tP') = 19.0GeV^/c'' and 
20.5GeVVc''. 

In the following, we illustrate the results of different 
fits using projected histograms of the slices of the Dalitz 
plot indicated by the vertical solid lines and horizontal 
dashed lines shown in Fig. [T] The three horizontal slices 
correspond to M^n'^'ip') regions below, around and above 
the Z(4430)+ mass region. The five vertical slices distin- 
guish the A'* (892) and M{KTr+) ~ I.4GeV/c2 regions 
and bands above, below and in between them. The sum 
of the latter three projections corresponds to the K* veto 



used in Ref. 0. 

FORMALISM OF THE DALITZ ANALYSIS 

The decay B A'tt+V' with the tp' reconstructed in 
the decay mode is described by four variables (as- 

suming the width of the ip' to be negligible). These are 
taken to be M{tt+iP'), M{KTr+), the ip' helicity angle 
(9) and the angle between the tp' production and decay 
planes ((/)). In this analysis we integrate over the angular 
variables 9 and cf>. The MC indicates that the reconstruc- 
tion efficiency is almost uniform over the full </> angular 
range; after integration over this angle the contribution 
from interference between the different tp' helicity states 
is negligibly small. This allows the ip' to be treated as a 
stable particle in the Dalitz analysis. 

In the V'' ^ TT+TT^J/-)/' channel, the tp' is likewise 
treated as stable. The tt+tt" system in this decay is 
predominantly in an S'-wave in this limit, the ■0' and 
J/ijj helicity states are the same, and we again find negli- 
gible interference contributions after integration over de- 
cay angles. Thus, our approach is the same as in the 
Dalitz analysis of the 5° K^-k^XcI decays in Ref. 0]. 

The amplitude for the three-body decay B — > Kir'^tp' 
is a sum over different quasi-two-body modes; resonances 
are described by relativistic Breit-Wigner functions with 
angular dependence. As the default fit model, we include 
all known low- lying Ktt'^ resonances [the k or AT* (800), 
and the K*{892), A:*(1410), A'^(1430), Ar|(I430), and 
Ar*(1680)] and a single exotic ir'^ip' resonance. In ad- 
dition to the physics model, the fit function includes a 
background term derived from the AE sidebands and 
is modulated by the MC-determined experimental effi- 
ciency. The MC sample is generated using the world- 
average -ip' branching fractions Q while to fix the relative 
fractions of the and B~^ contributions we use isospin 
symmetry. The Dalitz plots for the AE sideband and the 
MC sample are smoothed. The expression for the ampli- 
tudes, signal component of the fit function, and other 
details of the fitting procedure are the same as used in 
the analysis described in Ref. 0]. 

FIT RESULTS 

The eight projected Dalitz plot slices with fit results 
for the default model superimposed are shown in Fig. [21 
The Z(4430)+ signal is most clearly seen in the third ver- 
tical slice. The sum of the 1st, 3rd and 5th vertical slices 
(i.e. a Dalitz plot projection with the K* veto applied) 
is shown in Fig. [3l The n'^^p' resonance parameters de- 
termined from the fit are M = (4443+i[2) MeV/c^ and 
r = (107l43)MeV. The central values agree well with 
the parameters reported in Ref. 0, while the errors are 
somewhat larger. The statistical significance, calculated 
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FIG. 2: Dalitz plot projections for the slices defined in Fig.[T] 
(a)-(e) correspond to vertical slices fi'om left to right, (f)-(h) 
correspond to horizontal slices from bottom to top; in (f)-(h), 
plots including the full vertical scale are shown inset. The 
points with error bars represent data, the solid (dotted) his- 
tograms are the fit results for the default model that includes 
all low- lying K-k resonances and a single (without any) tt"'"?/)' 
state, and the dashed histograms represent the background. 
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M^(n*V'), GeV^/c" 

FIG. 3: The Dalitz plot projection with the K* veto applied. 
The points with error bars represent data, the solid (dotted) 
histogram is the Dalitz plot fit result for the fit model with 
all K-K resonances and a single (without any) -k^^' state, and 
the dashed histogram represents the background. 



from the change in 2 \ogC when the Z(4430)+ is included 
in the fit (taking the added degrees of freedom into ac- 
count) is 6.4 cr. The fit fractions and significances for all 
of the components are listed in Table HI The confidence 



TABLE I: The fit fractions and significances of all contribu- 
tions for the fit models with the default set of K-r'^ resonances 
and a single tt^V' resonance. 



Contribution 



Z(4430) + 

if* (892) 
K*{1A10) 
7^0* (1430) 
A-2*(1430) 
;^*(1680) 



Fit fraction (%) 



Significance 



5.7; 
4.1 



-1-3.1 
1.6 

+3.4 



64.81^*5 

5.3 ±2.6 
5.51* 
2.8"* 



- + 1.6 

-1.4 

h5.8 
-1.0 



6.4 a 

1.5 o- 
large 
0.5 a 
1.3a 
3.1a 
1.2a 



level (C.L.) of the fit model with (without) the Z(4430)+ 
is 36% (0.1%). The C.L.'s are determined using ensem- 
bles of the MC simulated experiments. 

To study the model dependence, we consider a variety 
of other fit hypotheses. These include: successively re- 
moving each K* resonance component; adding, for each 
case, a non-resonant phase-space term; relaxing the con- 
straints on the K mass and width; replacing the k with 
the LASS group's parameterization for the Kn S-wave 
amplitude [llj, and including another J = 1 ( J = 2) K* 
resonance with mass and width left as free parameters. 
The lowest ^(4430)+ significance of 5.4 ct corresponds to 
the model with a non-resonant phase-space term and a 
new J = 2 K* resonance. We treat the maximum varia- 
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tion of the Z{4430)~^ parameters from these different fit 
models as the systematic uncertainty. The resulting un- 
certainty estimates are given in the first row of Table [TTl 

We find the uncertainty due to the variation of the r 
parameter in the Blatt-Weisskopf form factors to be 
negligible. The contribution of the uncertainties in the 
mass and width of intermediate K* resonances that are 
fixed in the fit is also found to be negligible. 

We vary the assumption about the value of the B de- 
cay orbital angular momentum (L) for those cases where 
several possibilities exist. The resulting uncertainties are 
given in the second row in Table [TTl 

In the fits described above, the spin of the Z{M30)'^ is 
assumed to be zero. We find that the J = 1 assumption 
does not significantly improve the fit quality. The vari- 
ations in the Z(4430)+ parameters for the different spin 
assignments are considered as systematic uncertainties 
and are listed in the third row in Table [III 

We consider alternative smoothing procedures for AE 
sidebands and MC samples. The corresponding variation 
of the Z(4430)+ parameters are given in the fourth row 
in Table ni 

To obtain the total systematic uncertainties, the 
values given in Table arc added in quadrature. 
The resulting mass, width and fit fraction arc AI = 
iAAA3tlltll)McY/c\ T = (107;»^+74) MeV, / = 

7+3.1+9.4\(y 

(,o-'_i.6-2.7J Ro- 



table II: Systematic uncertainties in the Z(4430)^ mass, 
width and fit fraction due to various sources. 



Fit model 

L assignment 

Z spin assignment 

Smoothing procedure 



M, MeV/c^ r, MeV Fit fraction, % 



+ 14 


+56 


+3.6 


-13 


-52 


-2.7 


+8 


+44 


+2.0 


-0 


-0 


-0.0 


+9 


+8 


+8.4 


-0 


-0 


-0 


+4 


+ 17 


+0.5 


-3 


-23 


-0.2 



OTHER FITS 

In principle, more complex mass structures can be 
produced by reflections from higher Ktt^ partial waves. 
To examine this, we perform the Dalitz plot fit with a 
1780) resonance term added to the default model 
(see Fig. |4|). In this case, the Z(4430)+ signal persists 
with mass and width within 1 cr of their default model 
values and with a statistical significance of 4.7 ct. How- 
ever, the i<:|(1780) fit fractions — 6.8% and 6.6% for the 
Z(4430)+ and non-Z(4430)+ hypotheses, respectively — 
are very large for a resonance with a peak mass that 



is - 180MeV/c^ (~ 1.2r 



^1(1780), 



above the kinematic 



limit for B KTr~^ip' decays and for which only a small 
portion of the low-mass tail of the resonance is accessible. 
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FIG. 4: Dalitz plot projections as described for Fig. [J] The 
points with error bars represent data, the solid (dotted) his- 
tograms are the fit results for the model that includes all 
low- lying Kn resonances as well as the ifa (1780) and a single 
(without any) -K'^tp' state, and the dashed histograms repre- 
sent the background. 
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Moreover B if|(1780)V'' decay has an L > 2 centrifu- 
gal barrier and the (1780) — > Kit branching fraction is 
only (18.8 ± 1.0)% For these reasons, the 4.7 cr signif- 
icance estimate from this fit model is likely to be an un- 
derestimate of the real value. Studies of other B decays 
where the if|(1780) can contribute (e.g. B -> K 7r+, J /ip 
and B — > K*'K^ip') may provide further insight. The 
C.L. of the default model with an additional contribu- 
tion from the i^|(1780) is 58% (6%) for the Z(4430)+ 
(non-Z(443G)+) hypothesis. The significant Z(4430)+ 
contribution is concentrated in a small area of the Dalitz 
plot. 

If a second Z state is added to the fit, we find a mass 
M ~ 4.3 GeV/c^ and a width T - 0.2 GeV, with a sig- 
nificance of 3.9 cr. 

Angular distributions for ip' decays can be predicted 
based on the Dalitz plot fit results and therefore provide a 
useful cross-check (see Ref. 0] for details) . We find good 
agreement between data and predictions for various fit 
models. The statistics are not sufficient to discriminate 
between models in our approach. 

BRANCHING FRACTIONS 

To measure branching fractions we use only B^ — > 
K~tt^iIj' decays. The yields of these decays with the 
ijj' reconstructed in the £^£~ and tt+tt" J/-0 channels are 
found from fits to the AE distributions to be 1089 ± 34 
and 1166 ± 37, respectively. 

To determine the experimental efficiency, we used the 
phase-space MC events weighted according to the results 
of the Dalitz plot fit. The efficiencies are (19.2 ± 1.4)% 
and (8.2 ± 0.7)% for V' and V' ^ tt+tt" J/V' 

channels, respectively. The uncertainties include the de- 
pendence on the Dalitz plot model (0.1%); data and 
MC differences for track reconstruction (1% per track), 
and particle identification (4% for the K~tt~^ pair and 
4.2% for e+e-)] and MC statistics (0.6%). The un- 
certainties from different sources are added in quadra- 
ture. The efficiencies are corrected for the difference in 
lepton identification performance in data compared to 
MC, (-4.5 ± 4.2)%, determined from J/iP £+£- and 
e+e" — > e+e~£~*"^~ control samples. 

Using (656.7±8.9) x 10^/2 as the number of B°B° pairs 
and world-average values for the intermediate branch- 
ing fractions Q, we determine BlB'^ K~Tr~^'ip') = 
(5.68 ± 0.13 ± 0.42) x lO"''. This value is in agree- 
ment with BaBar's value (5.57 ± 0.16) x 10"'' 01 (sta- 
tistical error only). The branching fractions calculated 
for the two ip' decay channels arc in good agreement: 
(5.73 ±0.18) X 10-4 and (5.62 ±0.18) x lO"'' (statistical 
errors only). The systematic error includes contributions 
from the uncertainties in the efficiencies and the branch- 
ing fractions of the intermediate resonances. 

Based on the ^(4430)+ fit fraction we find a prod- 



uct branching fraction B{B° if-Z(4430)+) x 

B(Z(4430)+ Tr+^') = {3.2tlitli) x lO'^. This is 
in agreement with the previous Belle result 0] and con- 
sistent with the BaBar upper limit of 3.1 x 10~^ 

The dominant feature of the B — > Kir^ip' decay pro- 
cess is the B K*{892)ip' intermediate state. Us- 
ing the fit fraction from Table U and evaluating sys- 
tematic errors using the same procedure as for the 
^(4430)+ measurements, we determine the branching 
fraction B{B° if*(892)V) = (5.52t°;|^t°:i) x lO"* 
and the fraction of A'* (892)° mesons that are longitu- 
dinally polarized /l = (44.812:7^5:3)%- The branching 
fraction is somewhat below the world-average value of 
(7.2 ± 0.8) X 10~* [1]; the longitudinal polarization frac- 
tion agrees with the CLEOII result of 0.45 ± 0.11 ± 0.04 
and has better precision [l2|. 

CONCLUSIONS 

From a Dalitz plot analysis of i? ^ Kn'^ip' decays, 
we find a signal for Z(4430)''' tt+i/;' with a mass 
M = (UAStlltll) MeV/c2, width F = {lOltlttlt) McV 
and product branching fraction B{B^ A'^Z(4430)+) x 
S(Z(4430)+ ^ TT+V/) = (3.2lJJl^;^) x lO'^. The sta- 
tistical significance of this signal is 6.4 ct; the significance 
including systematic uncertainty from the fit models is 
5.4(7. These results agree with, and supersede previous 
measurements based on the same data sample reported 
in Ref. 0. 

In addition we determine the branching fraction 
B{B° -> if*(892)V) = (5-52ll]:ii!]:i) X 10-4 and the 
fraction of AT* (892) mesons that are longitudinally polar- 
ized Jl = (44.8^2:715:3)%- These are the first measure- 
ments of these quantities that are derived from a Dalitz 
plot analysis. 
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